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Exact solutions are obtained for the mean-field spherical model, with or 
without an external magnetic field, for any finite or infinite number TV 
of degrees of freedom, both in the microcanonical and in the canonical 
ensemble. The canonical result allows for an exact discussion of the loci 
of the Fisher zeros of the canonical partition function. The microcanonical 
entropy is found to be nonanalytic for arbitrary finite N. The mean-field 
spherical model of finite size N is shown to be equivalent to a mixed 
isovector/isotensor er-model on a lattice of two sites. Partial equivalence 
of statistical ensembles is observed for the mean-field spherical model in 
the thermodynamic limit. A discussion of the topology of certain state 
space submanifolds yields insights into the relation of these topological 
quantities to the thermodynamic behavior of the system in the presence 
of ensemble nonequi valence. 

Key Words: Phase transitions, spherical model, microcanonical, canon- 
ical, ensemble nonequi valence, partial equivalence, Fisher zeros of the 
partition function, RP W_1 cr-model, mixed isovector/isotensor cr-model, 
model equivalence, topological approach. 



1. INTRODUCTION 

Why discussing in detail an already simple model under even more sim- 
plifying assumptions? The spherical model is a simple model, constructed 
such as to be exactly solvable. In this article, this model is discussed under 
even more simplifying assumptions, namely mean-field type interactions, 
coupling each particle to each other at equal strength. The answer to the 
above question is manifold: 



1 Physikalisches Institut, Lehrstuhl fur Theoretische Physik I, Universitat Bayreuth, 
95440 Bayreuth, Germany; e-mail: Michael.Kastner@uni-bayreutli.de 

2 Institut fur Theoretische Physik III, Friedrich-Alexander-Universitat Erlangen- 
Nurnberg, Staudtstrafie 7, 91058 Erlangen, Germany 



1 



2 



Michael Kastner and Oliver Schnetz 



(i) The simplicity of the model allows for an exact solution, with or 
without external magnetic field, not only in the thermodynamic 
limit, but also for any finite number N of degrees of freedom, both 
in the microcanonical (section EJ) and in the canonical ensemble (sec- 
tionEJ). From these solutions, the calculation of finite-size quantities 
like the zeros of the canonical partition function is possible (section 
ED- 

(ii) With the exact result for the microcanonical entropy density sn, 
an example is at hand showing explicitly the nonanalyticity of sn 
as a function of the energy for arbitrary finite N (section EJ). This 
behavior is contrary to the well-known analyticity of the canonical 
free energy for finite systems. 

(iii) Although the microcanonical and the canonical ensemble describe 
different physical situations, both ensembles are known to yield 
equivalent results in the thermodynamic limit N —>■ oo under suit- 
able conditions (like the sufficiently short-rangedness of the inter- 
action potential).^ Systems with mean-field type interactions can 
be viewed as a caricature of (physical) long-range interactions, and 
they can serve as a playground for the investigation of ensemble 
nonequivalence.^ The mean-field spherical model is found to dis- 
play a peculiar type of nonequivalence, a phenomenon termed partial 
equivalence of ensembles (section |SJ). 

(iv) Two different statistical mechanical models are considered equiv- 
alent when their thermodynamic potentials are equivalent. The 
known examples of equivalent pairs of models are systems which, 
beside their differences in the Hamiltonian function, have important 
aspects in common, like lattice dimensionality and (infinite) system 
size. In section the mean-field spherical model for finite N is 
shown to be equivalent to a one-dimensional model, the so-called 
mixed isovector/isotensor cr-model on two lattice sites. This result 
implies equivalence of the mean-field spherical model with zero mag- 
netic field to the BP N ~ 1 cx-model. 

(v) In a series of p apers that appeared in the last 
years, ^ * ^^^ 1 ^ 11 * 12 * 1 ^ topology changes within a family of 
certain level sets of the state space of a statistical mechanical model 
were found to be intimately connected to phase transitions. In 
section El exact results on topological quantities of the mean-field 
spherical model are confronted with the thermodynamic results, 
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clarifying the relation between this topological approach to phase 
transitions and thermodynamics in the presence of ensemble 
nonequivalence . 3 

Altogether, its simplicity and the long-rangedness render the mean-field 
spherical model an instructive and profitable laboratory for the investiga- 
tion of a variety of effects in the field of the statistical mechanics of phase 
transitions. 

2. MEAN-FIELD SPHERICAL MODEL 

The spherical model of a ferromagnet was introduced by T. H. Berlin 
and M. Kac' 14 ' in 1952. It is devised such as to mimic some features of 
the Ising model, while, in contrast to the Ising model, having a continuous 
configuration space. The Hamiltonian function of this model is given by 

{A — > R, 
a ^ -- Jij o-iCTj - h di, 
* i,j=l i=l 

where a state a = (ai, . . . , <7jv) is determined by the values of the N degrees 
of freedom <jj G R, i = 1, . . . , N. The coupling matrix determines 
the strength of the interactions between the i-th and the j'-th degree of 
freedom, and h is a spatially homogeneous external magnetic field. The 
state space A is an (A^-l)-sphere with radius VN, and it is for this reason 
that the model is called "spherical" . 

Berlin and Kac considered the case of ferromagnetic nearest-neighbor 
interactions, where = J > for all degrees of freedom being nearest 
neighbors on a (^-dimensional cubic lattice, and = otherwise. They 
obtained exact expressions for the canonical free energy density / as a 
function of the temperature T. A phase transition takes place for d ^ 3, 
as f(T) is found to be a nonanalytic function, whereas no transition occurs 
for d — 1 and d = 2 (for a review of results on the spherical model see 
Ref. 03). 

Even in the presence of long-range interactions, decaying algebraically 
with some power r~ q (q > 0) of the distance r on the lattice, an exact 
analysis of the spherical model is possible in the thermodynamic limit. ^ 
The existence or absence of a phase transition, as well as the values of the 

3 Due to their relevance for the interpretation of the topological approach to phase 
transitions, part of the results discussed in section|H|has already been reported (without 
derivations) in Ref. 
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critical exponents in case of a transition, are found to depend on both, 
the spatial dimension and the value of q. 

We consider a mean field-like simplification of the spherical model 
where all degrees of freedom interact with each other at equal strength. 
The Hamiltonian function of this model can be written as 



A — > R, 

"2N 



H: { 1 (N y n (2.2) 



i=l J i=l 



The state space A is, again, an (AT-l)-sphere with radius \N, and the 
coupling constants J^- = 4 = 1,...,N) were chosen system-size 
dependent in order to guarantee the extensivity of the energy. 

In the following we will discuss this model within the microcanon- 
ical and the canonical ensemble. Note that the terms "microcanonical" 
and "canonical" have been used in connection with the spherical model 
in a different context,^ referring to the question whether the spherical 

N 

constraint ^ of = N, confining the states a onto the (N — l)-sphere A, 

i=l 

holds strictly (microcanonically), or only on average (canonically). More 
appropriately, these two cases should be regarded as different models in- 
stead of different statistical ensembles. In the present article, the notions 
microcanonical and canonical will be used in their meaning from ensem- 
ble theory, referring to the physical situations of fixed energy and fixed 
temperature, respectively. 



3. MICROCANONICAL SOLUTION 

The fundamental quantity for equilibrium statistical mechanics in 
the microcanonical ensemble is the microcanonical partition function or 
density of states 

Sl N (e) =Aj r 1 J---J do-x.-Ao-NSiHia) - Ne] (3.1) 



as a function of the energy density e and indexed by the number of 
degrees of freedom. Here, 5 denotes the Dirac distribution and 

A N = ■■■ da v ..da N (3.2) 
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is a normalization constant. The iV-fold integrals extend over the whole 
state space A. We have 



2 N 



(3.3) 



,i=i 



and hence, for the mean- field spherical model (|2.2|) . we can write 



Q N {e) = A~ N l 



dcri...dcr 



N 



x S 



where 
A N -- 



^ / N \ 2 w 

\i=l / t=l 



Ne 





' N 


2VN5 









(3.4) 



iV 



d ( T 1 ...da JV 2v / iV5 XX ~ N 



j=i 



27r N/2 N (N-l)/2 

r(iv/2) 



(3.5) 



is the surface of an (N — l)-sphere of radius \N, and T denotes the 
(complete) Gamma function. The second Dirac distribution in ()3.4|) is 
the spherical constraint which restricts the states a to the state space A. 
By the properties of the Dirac distribution, the integrand in ()3.4j) can be 
written as 



<3>(h 2 -2e) 
Vh 2 - 2e 



N 



N 



(3.6) 



j=l \i=l 

where O denotes the Heaviside step function and 



i=l 



h l/2 (e) = h±Vh 2 - 2s. 



(3.7) 



After interchanging the sum over j and the integral, we shift cr, — > 



hj(e) and obtain 



Q{h 2 -2e) 



X 



A N VW^2ej-t 

x 2\fN8 



d(7i...d(7jv-l 



'JV-l 



N-l 



i=l 



i=l 



(3i 
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by integrating over a^. The argument of the Dirac distribution contains 
an (N — 1) x (N — 1) bilinear form, 



N-l 



i=l 



i=l 





T 


( ai 


\ 




A 






\CTjv-1 j 




\<?N- 





(3.9) 



where A is an (N — 1) x (iV — 1) matrix with entries Au = 1 + Sk,e, and 5 
denotes Kronecker's delta. The eigenvalues of A are N and 1, where the 
latter is (N — 2)-fold degenerate. Hence, by an orthogonal transformation 
of variables in ()3.8|) and a rescaling in the direction of the eigenvalue N, 
we obtain 



Q(h 2 - 2e) 



+oo 



E 



A Ny /N(h? - 2e) 

x 



do"i...dcrjv-i 



'7V-1 



1=1 



(3.10) 



For /i^(e) < 1, the integral in ()3.10|) describes the surface of an (N — 2)- 



sphere of radius J N (l — (and vanishes otherwise). With equation 

f!3.5|) we obtain 

r(f) e(/i 2 -2 £ ) J 



^) 



. , JV-3 

"£[l-h%e)]^e[l-h%e)} (3.11) 



r 



N-l 



) Ay tt(/i 2 - 2e) 



as a final expression for the microcanonical partition function of the mean- 
field spherical model (|2.2|) . where /ij is defined in (|3.7j) . 
The microcanonical entropy density is defined as 



Its thermodynamic limit, 



s N (e 



s e 



AT 



In ^(e). 



lim sjv(e), 

N—*oo 



(3.12) 



(3.13) 



is the microcanonical counterpart of canonical thermodynamic potentials 
like the Gibbs free energy density, and it constitutes the starting point 
for microcanonical thermodynamics. Making use of Sterling's formula, we 
obtain from the microcanonical partition function (J3.11|) 



s e 



In 



1 - (\h\ - Vh 2 -2ey 



(3.14) 
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Figure 1: Domain of the microcanonical entropy density s(e) for the mean- 
field spherical model ([2,2)1 in dependence of the value of the external magnetic 
field h. Points within the shaded region obey —1 < \h\ — Vh 2 — 2e < 1 and 
correspond to accessible states. 



defined for all values of e and h for which the radical and the argument 
of the logarithm are positive. For visualization, in figure Q the domain of 
s is illustrated in the (e, /i)-plane. Exemplary plots of the graph of s with 
and without external magnetic field h, respectively, are shown in figure 121 



4. CANONICAL SOLUTION 

The fundamental quantity for equilibrium statistical mechanics in the 
canonical ensemble is the canonical partition function 

+oo 

Z N {(3) = A- N Y J ■■■ J Ao-L.AoNe-P 11 ^ = N J den N (s) e^ 6 (4.1) 

A —oo 

as a function of the inverse temperature (3 = > 0, where ks is Boitz- 
mann's constant and T is the temperature. The normalization constant 
An is the same as in ([3.2)1 . and the integration, as in the definition of 
the microcanonical partition function (|3.1|) . extends over the whole state 
space A. Inserting the microcanonical partition function (|3.11|) into ()4.1|) 
and interchanging the sum over j with the integral, the canonical partition 
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Figure 2: Exemplary plots of the graph of the micro canonical entropy density 
s for the cases h = and h = 1. The domain of s is the interval (— ^,0] for 

h = 0, and (— |, ^) for /i = 1. Note that, for h = 0, the slope 9s q^ ^ 1 has a 
positive lower bound, whereas the slope is unbounded for all h ^ 0. 



function reads 

Zn ^ = r( A 2 )v ^ J dz (1 - * 2 ) 2 ex P (! + *)} ( 4 -2) 

where z = ^-(e) has been substituted. This expression can be written 
as an infinite sum of confluent hypergeometric functions \Fi. For zero 
external magnetic field h = 0, a single summand survives, yielding 

Z h N =°(P) = 1 F 1 (if,2f) (4.3) 

for the zero-field canonical partition function of the mean-field spherical 
model. 

The canonical free energy density is defined as 



fs{fi) = ~^Z N {0). (4.4) 



Its thermodynamic limit, 



/(/?) = lira f N ((3), (4.5) 

N— >oo 



constitutes the starting point for canonical thermodynamics. An exact 
expression for f(/3) can be obtained in various ways: either by a large 
iV-expansion u& of the integral in (|4.3J1 . or, analogously to the calcu- 
lations of Berlin and Kac,' 14 ' by making use of the method of steepest 
descent. Having already at hand an expression for the microcanonical en- 
tropy density s, the free energy density is most easily obtained by means 
of a Legendre-Fenchel transform, 

f((3) = M[e-(3- 1 s(e)}. (4.6) 
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Evaluating this expression, we find that 



f{P) = e{p)-pr 1 8[E{p)] 



(4.7) 



with the following expressions for e((3): 

For the case of zero external magnetic field h = 0, 




(4.8) 



yielding 



for p ^ 1, 

for/?>l, 



(4.9) 



for the canonical free energy density for zero external magnetic field. The 
cusps in the graphs of e{(3) and f((3) can be viewed as a consequence of 
the fact that, for h — 0, the microcanonical entropy density (|3.14j) has a 
compact support and its slope has a positive lower bound (see figure |2J). 

For nonzero field h, e((3) can be expressed as the real root of a third 
order polynomial in e, 



8(3 2 e 3 - A(3{(3h 2 -2(3 + 2)e 2 - 2(<o(3 2 h 2 - (3 2 - 2(3h 2 + 2(3- l)e 

+ (3h 2 (4(3h 2 - (3 + 2) = 0. (4.10) 



In both cases, the expressions for e((3) are derived from the inverse function 
of the positive part [(3(e) > 0] of the microcanonical inverse temperature 



Exemplary plots of the graph of / are shown in figure 01 

Physically, the observed nonanalyticity in the zero field case marks 
a transition from a ferromagnetic phase with nonvanishing magnetization 
at low temperatures ((3 > 1) towards a paramagnetic phase with zero 
magnetization at high temperatures ((3 < 1). 

5. EQUIVALENCE TO AN ISOVECTOR/ISOTENSOR cr-MODEL 

Two different models are considered equivalent when they have iden- 
tical thermodynamic potentials. Well-known examples of systems being 
equivalent in the thermodynamic limit include 



ds(e) _ 1 + 2e - (1 - 2e)h(h 2 - 2s) 
de ~ (1 + 2ef - Ah 2 



1/2 
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Figure 3: The graph of the canonical free energy density /, showing a nonan- 
alyticity at (5 = 1 for h = 0, whereas / is smooth for h ^ 0. 



(i) a lattice gas of (quantum) particles, whose grand canonical poten- 
tial equals the canonical thermodynamic potential of an (anisotropic 

quantum) spin-| system in an external magnetic field on the same 
lattice,*™ 

(ii) the D- vector model (comprising the Ising, the XY-, and the Heisen- 
berg model as cases D = 1, 2, 3, respectively) which, on a hypercubic 
lattice and in the limit D — > oo, is equivalent to the spherical model 
on the same lattice. EEU 

Note, however, that in general such equivalence does not extend to other 
quantities apart from thermodynamic potentials and quantities derived 
therefrom.^ 

A more intricate relation is found for the mean-field spherical model 
and an extension of the RP D_1 cr-model, the so-called mixed isovec- 
tor/isotensor cr-model. The PJP D_1 a-model is a lattice model with an or- 
der parameter taking values in the real projective space HP 15-1 . This sys- 
tem has attracted some interest as a model of nematic liquid crystals,^ 
as well as in the particle physics community as a laboratory for investi- 
gating asymptotic freedom, topological phenomena, and universality. For 
a list of references and an account of the history of the PP D_1 a-model, 
the reader is invited to consult the article by Sokal and Starinets.^ In 
the same reference, the following extension of this model is discussed. 

The mixed isovector/isotensor a-model is defined by the lattice 
Hamiltonian 

((pP- 1 )" — > R, 
H : W ^-DZ [W ^ + | fa ■ Oi) 2 ] , ^ 

where the sum runs over all pairs (each pair counted once) of 

nearest neighbors on a lattice. Each component Cj of a state vector 
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a = (o"i, . . . , (Xjv) represents a point on a (D — l)-dimensional unit sphere 
§£>-i_ For = o, (jSHJ) reduces to the Hamiltonian of the RP ' 1 a- 
model. In reference , the canonical free energy density of the mixed 
iso vector /isotensor cr-model for arbitrary D on a one-dimensional lattice 
of size N with open boundary conditions was shown to be of the form 



J N,D 



r 



■ 2 ■ 



r 



D-l 



dz {i-z 2 y 



exp \ f3Dz [h + - 



N-l 



(5.2) 

Upon setting N = 2 in this expression and identifying D with N, the 
canonical partition function ()4.2|) of the mean-field spherical model is re- 
covered and the two models under consideration are found to be equiv- 
alent. Remarkably, and unlike the cases of model equivalence mentioned 
above, the models differ in fundamental aspects: The mean-field spherical 
model of finite size N with an external magnetic field is equivalent to the 
mixed isovector /isotensor cr-model with D = N on a (one-dimensional) 
lattice of only two sites. For the sake of clarity, the characteristics of the 
two equivalent models are confronted in the following table: 



model: 


mean-field spherical 


RP^" 1 a 


lattice dimension d: 
system size: 

degrees of freedom per site: 


N 
1 


1 

2 

N - 1 



Note that, for the two models considered above, model equivalence can be 
established even on the level of the Hamiltonian functions. The Hamil- 
tonian (|5.1|) of the mixed isovector /isotensor cr-model consisting of two 
lattice sites can be mapped onto the Hamiltonian (|2.2|) of the mean-field 
spherical model by first rescaling o\ and o"2 by a factor of D~ 1 ^ 2 , and then 
exploiting the rotational symmetry of the model by setting a\ = (!,...,!). 



6. FISHER ZEROS OF THE CANONICAL PARTITION FUNCTION 

The rare situation that a statistical system of N interacting degrees 
of freedom can be solved for arbitrary N allows to discuss the exact loca- 
tion of the so-called Fisher zeros, ^ the zeros of the canonical partition 
function in the complex temperature plane. The Fisher zeros form the 
natural analogon to the zeros of the grand canonical partition function in 
the complex fugacity plane as introduced by Yang and Lee in their seminal 
article. ^ 

For the mean-field spherical model, a discussion of the Fisher zeros of 
the zero-field canonical partition function (J4.3|) amounts to analyzing the 
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Figure 4: Szego curve (|6.ip in the complex /3-plane, together with some zeros 
of f , ^) for N = 1000 (circles) and N = 2000 (squares). 



complex zeros of the confluent hypergeometric function iFx(h, y, •) with 
respect to its third argument. Focusing on the large N asymptotics, we 
are interested in the complex zeros of y, j ) with respect to f3 in 

the limit of large N. At this point, we can benefit from the work by Sokal 
and Starinets:^ They proved that, as illustrated in figure the complex 
zeros lie close to the two branches with real part dt(/3) > 1 of the Szego 
curve 

| /3 |2 = e 2 W )-l) ) (61) 

approaching the curve and the value (3 — 1 in the thermodynamic limit 
of large N. This is a beautiful illustration of the Yang-Lee-mechanism of 
analyticity breaking, in which, for large but finite N, the complex zeros of 
a real-analytic function (the canonical partition function) approach and 
pinch the real axis, giving rise to a nonanalyticity in the limit N —>■ oo. 
More precisely, in this limit the zeros lie dense on the branches of the Szego 
curve, thus cutting the complex plane into two disconnected domains. 

7. NONANALYTICITY OF THE FINITE-SYSTEM ENTROPY 

As mentioned in the preceding section, the canonical partition func- 
tion, and hence also the canonical free energy density as defined in 
(|4.4jl . are known to be real-analytic functions of (3 for any finite number 
iV of degrees of freedom. For microcanonical thermodynamic functions 
like the entropy density sn, an argument guaranteeing analyticity does 
not exist. Nonetheless, the belief that sn were real-analytic in e, serving 



On the mean-field spherical model 



13 



as a justification for certain arguments or operations, has been expressed 
in several publications.^^ 

With the result 1)3.11)1 for the microcanonical partition function Qn of 
the mean-field spherical model, we have an exact result at hand from which 
it is straightforward to read off the nonanalyticity of Qn, and therefore of 
sn, for all values of the external magnetic field h satisfying < \h\ < 1. 
The nonanalyticity is located at e = \h\ — ^, as at this value the number 
of nonvanishing contributions in the sum over j in (j3.11)) jumps from one 
to two. The nonanalyticity is such that a discontinuity occurs in the A-th 
derivative of sn(e) with respect to e, where 

(7.1) 

i.e., the largest integer smaller than or equal to — jp. The observation 
that A is of order N for N — > oo is in agreement with ()3.14)1 where we 
found a smooth microcanonical entropy density s in the thermodynamic 
limit. 

One might argue that such a nonanalytic behavior of s is of little 
physical relevance for the system sizes N under consideration in statistical 
physics. For suitable choices of the interaction potential, however, we 
expect to find models of physical interest for which A is of order 1.^3 

8. NON EQUIVALENCE OF STATISTICAL ENSEMBLES 

As already pointed out in the introduction, the sufficiently short- 
rangedness of the interaction potential is a sufficient condition for the 
equivalence of the statistical ensembles in the thermodynamic limit. If, 
however, the microcanonical entropy density s is explicitely known, more 
can be said: the equivalence, or nonequivalence, of statistical ensembles 
can, loosely speaking, be traced back to the concavity properties of s (see 
the article by Ellis, Haven, and Turkington^ for a precise formulation of 
this statement as well as of the definition of ensemble equivalence). 

For a third type of relation of the ensembles, situated somewhere 
in-between equivalence and nonequivalence, the term "partial equiva- 
lence" has been coined,^ and it is this kind of situation we encounter 
for the mean-field spherical model with zero external magnetic field 
h = 0: The microcanonical entropy density, as depicted in figure |2J is 
a strictly concave function with compact support. Furthermore, we find 
f3{e) = G [l,oo), i.e., a positive lower bound exists on the micro- 
canonical inverse temperature. As a consequence of these properties and 
of the fact that the canonical free energy density / can be derived from 



A = X(N) 



N-2 
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s be means of a Legendre-Fenchel-transform ()4.6|) . the following relation 
arises: For all values of (3 ^ 1 above the lower bound on the microcanonical 
inverse temperature, a one-to-one correspondence between microcanonical 
and canonical macrostates exists, and the two ensembles are said to be 
equivalent. For (3 < 1, however, all the various canonical macrostates cor- 
respond to the microcanonical macrostate with maximum energy e = 0. 
This is clearly not a one-to-one correspondence and the ensembles are 
nonequivalent for this range of inverse temperatures. It is in this sense 
that, in spite of the strict concavity of the microcanonical entropy den- 
sity s, the equivalence between the ensembles is only "partial" for the 
mean-field spherical model. 

For arbitrary nonzero values of the external magnetic field h ^ 0, the 
microcanonical entropy density is strictly concave and is unbounded 
below and above, hence equivalence of the microcanonical and the canon- 
ical ensemble holds. 

Focusing on the question whether a phase transition is present in 
one ensemble or the other, the following considerations are of interest: 
one of the common choices is to define a phase transition as the presence 
of a nonanalyticity in the canonical free energy density. For all finite 
N, this quantity is known to be a real-analytic function of the (inverse) 
temperature. Only in the thermodynamic limit iV — ► oo a nonanalyticity, 
and hence a phase transition, may be found. This is exactly the situation 
we observe for the mean-field spherical model with zero external magnetic 
field h = 0: For systems consisting of a finite number N of degrees of 
freedom, the canonical free energy density /jv(/3) is the logarithm of a 
confluent hypergeometric function with all positive arguments [equations 
(14. 3|) and (J4.4|) ]. and is therefore real-analytic. In the thermodynamic 
limit, however, the nonanalyticity in f(f3) at (5 = 1 is obvious from (J4.9)) . 

In the microcanonical ensemble, conditions on the microcanonical 
entropy density s for the existence of a phase transition are somewhat less 
established. An obvious way to start is to look for conditions on s which, 
in the case of ensemble equivalence, correspond to the definition of a phase 
transition from the canonical free energy density. 4 Such conditions on s(e) 
are 

(i) the existence of an open interval of energies e for which s(e) is not 
strictly concave, or 

(ii) a nonanalyticity in s(e). 

4 A different definition of phase transitions, based on the microcanonical entropy, is 
proposed in Ref. ^ . 
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Then it might appear reasonable to adopt the same conditions also for 
systems for which equivalence of ensembles does not hold. Interpreting the 
microcanonical entropy density s obtained for h = in the thermodynamic 
limit in 1)3.14)1 along the lines of (i) and (ii), it is straightforward to discuss 
the issue of existence of a phase transition: s is a real-analytic and strictly 
concave function on its entire domain, therefore no phase transition takes 
place in the microcanonical ensemble. 

The situation observed for the mean-field spherical model with h = 
is typical for systems showing partial equivalence of ensembles. (For fur- 
ther examples see reference ^.) In spite of the fact that partial equiva- 
lence appears to be a relatively "mild" form of ensemble nonequivalence, 
its impact on the existence of phase transitions can be substantial, result- 
ing in the existence of a phase transition in the canonical ensemble, and 
the absence in the microcanonical one. 



An entirely different approach to phase transitions, not making use 
of any of the Gibbs ensembles, has been proposed a couple of years ago. 
This topological approach connects the occurrence of a phase transition to 
certain properties of the Hamiltonian function H, resorting to topological 
concepts. From a conceptual point of view, this approach has a remarkable 
property: The microscopic Hamiltonian dynamics can be linked via the 
maximum Lyapunov exponent to the topological quantities considered.^ 
With the topological approach, in turn, linking a change of the topology 
to the occurrence of a phase transition, a concept is established which 
provides a connection between a phase transition in a system and its 
underlying microscopic dynamics. 

The topological approach is based on the hypothesis'^ that phase 
transitions are related to topology changes of submanifolds S e of the state 
space of the system. The S e consist of all points q of the state space A 
for which the energy H(a) per degree of freedom equals a certain level e, 
i.e., 



(Or, in a related version, the topology of submanifolds M £ consisting of all 
points q with H(q)/N ^ e is considered.) The hypothesis then conjectures 
that a topology change within the family {S e } at e = e c is a necessary 
condition for a thermodynamic phase transition to take place at e c [or at 
the corresponding critical temperature T c = T(e c )]. This hypothesis has 
been corroborated by numerical and by exact results for a model showing 
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a first-order phase transitional^* as well as for systems with second-order 
phase transitions. < 4 l 9 l 3:j l^t A major achievement in the field is the recent 
proof of a theorem, stating, loosely speaking, that, for systems described 
by smooth, finite-range, and confining potentials, a topology change of 
the submanifolds S £ is a necessary criterion for a phase transition. 
Recent results on the state space topology of the mean-field 4 -model 
have cast some doubt on the general relevance of topology changes for the 
occurrence of a phase transition, t^iMI Dl rt these findings require further 
investigation. 

A topology change in S £ is clearly not sufficient to entail a phase 
transition. This follows for example from the analytical computation of 
topological invariants in the fc-trigonometric model, ®^ where the number 
of topology changes occurring is shown to diverge in the limit N — > oo, 
although only a single phase transition takes place. So topology changes 
appear to be rather common, and only particular ones are related to phase 
transitions. While some ideas relating the "strength" of a topology change 
to the occurrence of a phase transitions have been put forward,^ there are 
strong indications that such a criterion based exclusively on topological 
quantities cannot exist in general.^ 

Ribeiro Teixeira and Stariolo have analyzed the topological struc- 
ture of the state space submanifolds of the mean-field spherical model in 
reference Completing their findings and translating them into our 
notation, the topology of the submanifolds S e can be written as 

( S N ~ 2 for -\h\ - \ < e < \h\ - §, 

S £ ~ \ S N - 2 + S N - 2 for \h\ - | < e < f and \h\ < 1, (9.2) 
y otherwise, 

where ~ indicates topological equivalence, the +- sign denotes a topological 
(unconnected) sum, and is the empty set. The regions of the (e, /i)-plane 
corresponding to the different topologies and the lines at which topology 
changes occur are plotted in figure |SJ The only topology change in the 
interior of the physical parameter domain is from S^ -2 to S^ -2 + S^ -2 . 
This change occurs for < \h\ < 1 and e = \h\ — | and is responsible for 
the nonanalyticity of s^(e) (see section Q). 

For finite N, we consider this connection between topology and ana- 
lyticity to be generic. It is a deeper question whether or not the nonana- 
lyticity survives the thermodynamic limit. The topological approach is an 
attempt to answer this question and to establish a connection between the 
topology of the state space submanifolds S e and thermodynamic quanti- 
ties of the system. For systems like the mean-field spherical model, where 
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Figure 5: Lines in the (e, /i)-plane where topology changes occur. The shaded 
region marks the support of the density of states f2jv(£), i.e., the values of e 
which, for a given external field h, are accessible for the system. The various 
regions in the plane separated by the bold lines correspond to the following 
submanifold topologies: (i) to the left of the shaded region: T, £ ~ 0, (ii) the 
two shaded triangles: S £ ~ S N ~ 2 , (iii) the small, triangle shaped shaded region: 
S e ~ S N ~ 2 + § N ~ 2 , (iv) to the right of the shaded region: £ e ~ 0. 



nonequivalence of the statistical ensembles is observed, one would expect 
a correspondence between topological and thermodynamic quantities only 
with regard to one of the statistical ensembles. 

Comparing the topology of S e in (J9.2J) with the canonical thermody- 
namics of the mean- field spherical model with zero external field h = 0, 
one would associate the phase transition observed in (|4.8|) at inverse tem- 
perature (3 = 1 and internal energy e = with the topology change from 
g7v-2 _|_ £ a j. ^ ^ _ However, the same type of topology 

change is present for all \h\ < 1 at (e, h) = (^-,h), whereas canonical 
thermodynamics asserts the absence of a phase transition for all h ^ 0. 
Hence, one would conclude that the information contained in the topol- 
ogy of E £ (or, in a related description, of M £ ) might not be sufficient as 
to distinguish between the presence or absence of a phase transition. 

However, a meaningful connection between the topology of S £ and 
the thermodynamics of the model can be established in the microcanonical 
ensemble. Within this framework, following the discussion of section |H1 
the zero field case is not special, as a phase transition occurs neither with 
nor without an external field. This is in accordance with the topological 
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results where the topology change is the same for all \h\ < l. 5 

The mean-field spherical model is not the only example for which 
the observed topology changes find a meaningful interpretation only when 
compared to microcanonical quantities. For the mean-field Xy-model 
as considered in references O anc l J3 ; the observed discontinuity in the 
modulus of the Euler characteristic is interpreted as being related to the 
phase transition of the model at zero external field h = in the canon- 
ical ensemble. However, similarly to the mean-field spherical model, the 
same topology change is present for h 7^ 0, although a phase transition 
does not occur in this case. Again, a meaningful interpretation of these 
findings should be possible when comparing topology to microcanonical 
thermodynamics. 

In addition to the above observations based on examples, a general 
remark on the relation between the topological approach and the statistical 
ensembles can be made. The microcanonical ensemble and the topological 
approach, different as they are, have the same quantity as a starting point: 
The microcanonical partition function (J3.1|) of a system, which is at the 
basis of the microcanonical ensemble, can be defined as the volume of the 
level sets (|9.1|) . Considering the topology of the very same level sets, we 
arrive at the fundamental quantity underlying the topological approach, 
and, in the light of this observation, the close connection between these 
two concept appears plausible. 

10. SUMMARY AND CONCLUSIONS 

The mean-field spherical model has proved to be an instructive and 
profitable laboratory for the investigation of a variety of effects and the- 
ories in the field of the statistical mechanics of phase transitions. Due 
to its simplicity, an exact solution can be obtained for any finite or infi- 
nite number N of degrees of freedom, even in the presence of an external 
magnetic field. 

An exact expression for the microcanonical entropy density s as a 
function of the energy density e was obtained. For external magnetic field 
h with < \h\ < 1, s is found to be nonanalytic at e — \h\ — | for arbitrary 
finite N. This is in contrast to the canonical case, where the free energy 
density fjq as a function of the inverse temperature (3 is real-analytic for 
all finite N. 

Note, however, that for a class of systems with non-confining potentials the infor- 
mation from the submanifold topology was found to be insufficient as to distinguish 
between the occurrence and the absence of a phase transition, regardless of the statis- 
tical ensemble considered.^ 
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Also for the canonical partition function Z N as a function of (3 an 
exact expression was given, allowing a discussion of the loci of the zeros of 
the partition function in the complex /3-plane. For h = and in the limit 
of large N, these zeros are found to lie dense on certain branches of the 
Szego curve, approaching the real axis at (5 — 1. The observed pinching 
of the real axis is a beautiful illustration of the Yang-Lee-mechanism of 
analyticity breaking in the thermodynamic limit. 

Equivalence of the mean-field spherical model of finite system size 
N to a mixed isovector/isotensor c-model on a lattice of two sites was 
established. This result implies equivalence of the mean-field spherical 
model with zero magnetic field to the RP^ -1 cr-model. 

By comparing, for h — 0, the microcanonical and the canonical 
solutions of the mean-field spherical model in the thermodynamic limit 
N — > oo, ensemble equivalence is found to be violated. Instead, a phe- 
nomenon termed partial equivalence is observed. This can be viewed as 
a consequence of the fact that the microcanonical entropy s is a strictly 
concave function with compact support, while its derivative, the micro- 
canonical inverse temperature (3(e) = has a positive lower bound. 

A complete characterization of the topology of the state space sub- 
manifolds S e of constant energy e is possible for the mean-field spherical 
model. The result allows for a discussion of the relation of these topologi- 
cal quantities to the thermodynamic behavior of the system. Based on this 
example, the connection between topological and thermodynamic quanti- 
ties in the presence of ensemble nonequivalence was discussed, observing a 
close relationship between the topological approach and the microcanon- 
ical ensemble. 
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